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A method for the nonperturbative calculation of scattering amplitudes and cross 
sections is discussed in the context of light-cone quantization. The Lanczos-based 
recursion method of Haydock is suggested for the computation of matrix elements 
of the resolvent for the light-cone Hamiltonian, from which the T-matrix can be 
constructed. The scattering of composite particles is handled by a generalization 
of a formulation given by Wick. 



1 Introduction 

Recently there has been a significant amount of progress in the nonperturbatb/£ 
calculation of bound-state properties in light-cone quantized field theories. Era 
Comparable progress in the calculation oLscattering states has been lacking, 
although techniques have been proposed. The ability to compute outcomes 
of scattering processes is not only of interest in its own right, but also usefiil 
for construction of renormalization conditions in bound-state calculations. □ 
Here we consijier a calculational scheme based on numerical approximations 
to resolvents □ and on a generalization of a Hamiltonian formulation for the 



scattering of composites. 

The construction ia^dpne in terms of light-cone coordinates, where t -\- z 
takes the role of time. QEl This is not a necessary choice, but it does bring 
several advantages and makes immecliate contact with work that uses these 
coordinates for bound-state problems. EJd The advantages include a simple vac- 
uum, well defined Fock-state expansions, and kinematical boosts. This last 
advantage is particularly important for the scattering of composites because 
one can compute the state of a composite particle in one frame and boost to 
other frames without again referencing the interaction. 
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The remainder of this paper is structured as follows: Section ^ contains a 
description of the notation used for light-cone coordinates and a formulation 
of two-body cross sections in terms of these coordinates. The key quantity in 
this formulation is a matrix element of a resolvent; a numerical method for its 
computation is discussed in Sec. |^ and applied to a nonrelativistic example in 
Sec. 0. An extension to the scattering of composites is given in Sec. ||. Finally, 
Sec. ^ provides a summary and indications of future work. 

2 Cross Sections in Light-Cone Coordinates 

We use the following notation for light-cone coordinates: 

=i±z, x± = {x,y), (1) 

where is taken as the time variable. Momentum components are similarly 
defined 

P^=E±p^, P± = iPx,Py). (2) 

A dot product is then written as 

p-x^^{p^x^ +p^x'^)-p±-x±, (3) 

so that p~ is conjugate to x~^ and is recognized as the light-cone energy. 
The other momentum components form the light-cone analog of the three- 
momentum 

P=(P+,P±)- (4) 

The evolution of a state is determined by the Hamiltonian V~', however, the 
term "light-cone Hamiltonian" is also usedH for the combination 

Hlc = P^V~ , (5) 
where is the total momentum. The fundamental eigenvalue problem is 

/i,,r2 I 2 

H^c\p) ^ -i^\p) . m^P\p)- (6) 

For two-body scattering A + B ^ C + D, the center-of-mass cross section 
can be written B 

where Mfi is the invariant amplitude extracted from the S matrix 

5/, = -f (2^)4^(4) „ ^ 
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with states normalized according to 

The flux factors are determined by the current matrix element 

{A\AA) = \{A\{J+ + J-)\A) =2Ea. (10) 

with Ea the energy of particle A. For fermions, the spinors are normalized by 
uu = 2m. From the optical theorem we also have 

2ImMii ^^^^ 



In the light-cone interaction picture the invariant amplitude can be 
related to a matrix element of a T operator. The S matrix can be written as 

Sf. = {f\U{x+,x+)\i) = im - 2mTf,SipJ ~ p-) , (12) 

with the T operator defined by 

T=V + V—-^——V, (13) 

p + le — V 

where 

V- = + V. (14) 
On removal of the momentum conserving delta function 

Tf,^tf,5{p^-p.) , (15) 

we obtain 

= ^^^^ 

We then need a means to compute tfi. In particular, we need to handle matrix 
elements of the second term in T which involve the resolvent (p^ +ie — V~)~^. 

3 Recursion Method for Resolvents 

We divide the matrix element of the T operator into two pieces, the Born term 
and the resolvent term 

Tf^^V^^+Tf^, r)f ^(/IV^— (17) 



To the latter we apply the recursion method developed by Haydocki in the 
context of condensed matter physics. This method is based on the Lanczos 
diagonalization algorithm Ell which converts the original Hamiltonian operator 
into an approximate tridiagonal representation. The tridiagonal form is easily 
inverted to yield the resolvent. 

Define an initial vector Iri*) — cV\i), with c a normalization constant, and 
carry out Lanczos recursionsEj 



bi = 0, a„ = {Un\V \Un) , {Un\u„) = 1 
&„+l|w„+l) = V^\Un+l) - an\Un) - 6„|w„-l) . 



(18) 



Although this could in principle be done analytically, the usual approach is to 
begin with a large but sparse matrix representation for in a, Fock basis. 
A standard way of constructing such a representation is discrete light-cone 
quantization (DLCQ). Ertl The vectors provide an orthonormal basis in 
which has the tridiagonal representation 



/ai &2 
&2 02 63 
63 03 64 

64 ■ • 
0.. 

V 



(19) 



and the forward matrix element of the resolvent is approximated by the upper 
left matrix element of the inverse 



(2) 



(20) 



Non-forward matrix elements can be obtained by using the block Lanczos pro- 
cedure in place of the ordinary one described above. E3 

Unfortunately, the discrete nature of the calculation does not allow this 

finite matrix representation to be useful in all respects. In the limit of e — > 

(2) 

0, the imaginary part of this approximation to will contain only delta 
functions at the discrete eigenvalues of the finite matrix approximation A. To 
properly model the contribution from the continuum, an extension must be 
made. This can be done in terms of a continued fraction representation 
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p, +ie- ai 



(21) 



Pi +i€ — a2- 
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by replacing the last level with a terminator, a model for an infinite fraction. lil 
A model can be generated in various ways; however, a logical choice in the 
present situation is to use the second Born approximation as a generator. The 
matrix element of the resolvent for the free Hamiltonian 

{^V^^--V\^) (22) 

is computed exactly and then in continued fraction form by iterating with 
Vq an equal number of times. The exact form dictates the structure of the 
remaining infinite fraction, which is used as a terminator for the full resolvent. 

4 Nonrelativistic Example 

As an example of how the terminated recursion method works, we consider 
nonrelativistic S-wave scattering from an exponential potential. This is a com- 
mon choice for testing numerical methods for scattering problems because it 
has an exact solution. 

The partial wave expansion for the T-matrix is 

{k'\TiE)\k) = -i^ ^(2; + l)r,(fc', fc; E)PiCk' ■ k) , (23) 



where 



Tiik',k;E)^ j drukn{r)V{r)uu{r) (24) 

+ / drukn{r)V{r)——^-—-V{r)uki{r) (25) 
J E + te- Hi 



and 

2m dr"^ 2mr 

The partial wave contribution T; is related to the phase shift Si by 



Ukiir) = rjiikr) , Hi = - — — + + Vir) . (26) 



Ti{k,k-e/2m) = -TT^e'^' sin (5/ or t&n5i = . (27) 

For S-wave scattering^of a particle of mass m from V{r) = -e"''/"/ ma 
the exact solution is lJ 

,5o(A:) = Arg(J2^fca(2%/2)r(l -I- 2ka)) - 2kalnV2. (28) 
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To solve this problem numerically, |j we introduce a discrete approximation 
with M interior points between r = and a radial cutoff R. The step size is 

- TT-^ ■ (29) 

Integrals are approximated by the trapezoidal rule, and derivative operators 
by finite differences. The radial cutoff is kept larger than 27r/fc in order to 
include one complete wavelength at the chosen energy. The terminator for the 
recursion is generated from the second Born approximation 

(2) _ f j,.^ '^^//'^„(0) 



T/'^ ~ / drdr'uk'iir')Vir')g^^>{r',r)Vir)ukiir) , (30) 



where 



(0) , miK , . , 

9Eiyr,r) = — —rrji{Kr<)hJ{Kr>) (31) 



and K = \/2mE. 

Figure |] displays the results of calculations done for k — 0.25/a, R = 
30a > 27r/fc ~ 25a, and M = 49, 99, 199, 399, and 799. They show good 
convergence properties. More importantly, they show that the imaginary part 
of the T-matrix element is calculated reliably. 



5 Scattering of Composites 

The formulation in Sec. || assumes implicitly that the particles involved in 
the scattering process are elementary. In a field-theoretic setting this is not 
typically the case for a nonperturbative calculation. The particles A, B, C, 
and D can be themselves complicated bound states of the full Hamiltonian. It 
is then difficult to determine the nature of the interaction V that appears in 
the T operator. 

To construct the format the T operator in this case, we generalize a 
construction given by WickQ to compute the S matrix 

Sab^cd = {CD - \AB+) (32) 

from scattering eigenstates of the light-cone Hamiltonian iJLC- Let and -B^ 
create one- particle solutions from the Fock vacuum |0), so that 

|A)=At|0), i7Lc|A) = ^^%t|k|A). (33) 

■^The method presented here is not by any means the best way to solve this particular 
problem. It is introduced purely for the purpose of testing the recursion method. 
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Figure 1: Phase shift for S-wave scattering from an exponential potential, as eoniputed by 
the reeursion method with termination based on the second Born approximation. Ar is the 
spacing used in trapezoidal approximations to the integral expression for the partial-wave T 
matrix element. 



The two-particle scattering eigenstates of Hj,c are then 

\AB±) = B^A) + ^ „ Vb\A) , (34) 

Sab ±ie- Hj^c 

where 

Vs^[H^c,B^]-^^^^±^B^ (35) 

and 

is the invariant mass of the scattering state and the Mandelstam s of the 
scattering process. The definition of Vb is chosen to ensure that 

Hi^c\AB±) = sab\AB±) . (37) 

The T operator is constructed by reducing the S matrix to 

SaB^CD = ScD,AB - 2-Ki5{SAB — SCD)Tl,Cfi ■ (38) 

We find 

TLCfi = P+Tfi = {CD - \Vb\A) , (39) 
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which can be expanded as 



P+Tf, = {C\VI ^ Vb\A) + {C\DVb\A) , (40) 

SAB + «e - i/LC 



with 



{C\DVb\A) = (C\Vb -r—. ^^T-^ V^\A) 

- ^TJ^ ^^TJ^ - ^LC 

+ {C\[D,Vb]\A). (41) 

The first term in ( ^ ) can be interpreted as a direct scattering of B by A. The 
first term in ( |4l] ) is the crossed contribution, where D is produced before B is 
absorbed. The second term in (|4l|) is-a generalized seagull contribution, which 



did not appear in Wick's formalism. 

To use these expressions efficiently, {C\DVb \A) should be computed explic- 
itly rather than work with the decomposition into crossed and seagull terms. 
This is because the decomposition introduces an inverse of the Hamiltonian. 
The inverse is unavoidable for the direct term, and here one would apply the 
recursion method discussed in the previous sections. For both terms one needs 
the action of Vb on the one-particle state \A). 



6 Concluding Remarks 

A method for the nonperturbative calculation of cross sections has been pre- 
sented in the context of light-cone quantized field theories. The key expression 
is (^0[), where the light-cone T matrix is given in terms of known operators 
and comgosite-particle eigenstates. The terminated recursion technique of 
HaydockB is to be used to handle the invetse of the light-cone Hamiltonian. 
Work on an application to a simple model is in progress. One direction to 
consider for possible improvements to the method is to explore variaticuial for- 
mulations analogous to the Schwinger and Kohn variational principlesll3 used 
in nonrelativistic scattering. 
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